Energy relaxation in disordered charge and spin density waves 
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We investigate collective effects in the strong pinning model of disordered charge and spin density 
waves (CDWs and SDWs) in connection with heat relaxation experiments. We discuss the classical 
and quantum limits that contribute to two distinct contribution to the specific heat (a Ct, ~ 
contribution and a ~ T" contribution respectively), with two different types of disorder (strong 
pinning versus substitutional impurities). From the calculation of the two level system energy 
splitting distribution in the classical limit we find no slow relaxation in the commensurate case 
and a broad spectrum of relaxation times in the incommensurate case. In the commensurate case 
quantum effects restore a non vanishing energy relaxation, and generate stronger disorder effects 
in incommensurate systems. For substitutional disorder we obtain Friedel oscillations of bound 
states close to the Fermi energy. With negligible interchain couplings this explains the power-law 
specific heat d, ~ observed in experiments on CDWs and SDWs combined to the power-law 
susceptibility x{T) ~ r"^+" observed in the CDW o-TaSa. 

PACS numbers: 71.45.Lr, 05.70.Ln, 63.50.-fx, 75.30.Fv 
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I. INTRODUCTION 



There exist many examples of systems showing 
slow relaxation and ageing: spin glassesi*^, disordered 
dielectric g'^i^i^ , supercooled liquids?., etc ... Charge den- 
sity waves (CDWs) and spin density waves (SDWs)i2i2ii£ 
show "interrupted ageing"—, meaning that there exists 
an upper bound Tmax to the relaxation times. The pro- 
tocol of ageing experiments in CDWs and SDWs is the 
following: the system is equilibrated at the temperature 
T (one waits for a time longer than Tmax)- At time t = 
the temperature is changed from T to T + AT where 
AT > is very small compared to T. The tempera- 
ture is kept constant until the waiting time t^, where it 
is brought back to T. The heat flows between the CDW 
or SDW sample and the cold source are recorded as a 
function of time. Ageing in the thermal response takes 
place even for very small values of AT. Since very small 
temperature variations are applied in the experiment it 
is reasonable to suppose that the size of the correlated 
objects does not evolve in time and that the thermal re- 
sponse is due solely to the variation in the population 
of metastable states. This can be contrasted with the 
coarsening dynamics where the size of correlated domains 
increases with time. 

In a recent worlii^ we applied the idea of dynam- 
ical renormalization groupi^ii^ to calculate the spec- 
trum of relaxation times of a model of disordered CDW 

or SDWi5ii6,i7,i8,i9,20,2i,22,23,24^ including interactions 

among bisolitons. A drawback of this approach^^ is that 
we supposed a coarsening dynamics following a quench 



from high temperature, a situation that is not realized in 
experiments, and we were not able to address the waiting 
time dependence of the relaxation time spectra. One goal 
of the present article is to address these issues that were 
left open in our previous work^^, and to put on a mi- 
croscopic basis the random energ y-li ke (REM-like) trap 
model that was proposed in Ref. \12L and inspired from 
trap models developed for glasses and spin glassesiiiS^. 

More specifically we show here that heat relaxation ex- 
periments can be described by assuming two types of de- 
fects (strong pinning and substitutional impurities), cor- 
responding to the "classical" limit where the CDW or 
SDW is viewed as a classical elastic medium with bisoli- 
tons generated by strong pinning impurities distributed 
at random, and to the "quantum" limit where solitons 
due to substitutional disorder interact quantum mechan- 
ically by excitations of the gaped background. The exis- 
tence of two effects is in agreement with the experimental 
observation that the low temperature out-of-equilibrium 
specific heat can be decomposed into three contributions: 
(i) the Cv ^ tail of a Schottky anomaly at very low 
temperature (typically for T < 100 300 mK; the upper 
bound depends on the amplitude of the 1/T^ contribu- 
tion); (ii) a Cv ^ T" power-law specific heat with a ~ 
0.3 1.2 at intermediate temperatures (0.1 < T < 1 K); 
and (iii) the "trivial" contribution of phonons ^ T^ 
at high temperature (T > 1 K). By Schottky anomaly 
we mean that the equilibrium specific heat of a two-level 
system with energies Eq and Eo + AE is given by 



(AE)' 



4T2 cosh^ (A£:/2r) ' 



(1) 



*E-mail: melin@grenobIe.cnrs.fr 

tU.P.R. 5001 du CNRS, Laboratoire convcntionne avec I'Universite 
Joseph Fourier 



having a maximum (called a Schottky anomaly) at 
Tmax — 0.416Ai?. The specific heat is approximately 
equal to Cv ~ (AE)'^ /AT'^ in the large temperature tail. 

Following Ref. [23, the contribution (i) is interpreted in 
terms of two-level systems due to strong pinning impuri- 
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ties. The contribution (ii) is interpreted as midgap states 
interacting through Friedel oscillations. Friedel oscilla- 
tions of a single impurity were probed directly by x-ray 
diffraction experiments in Ref. \2&. Another evidence in 
favor of the coexistence of strong pinning and substitu- 
tional impurities is that the CDW compound o-TaSs can 
be doped by Nb, a substitutional impurity. This changes 
only the amplitude of the Cy ^ T" contribution, but 
leaves unchanged the Cy ~ contribution^i, suggest- 
ing that the power-law contribution is related to substi- 
tutional disorder. Even though not affected by substi- 
tutional disorder we do not interpret the ~ l/T^ as 
a property of the pure compound. A nuclear hyperfine 
interaction can be excluded from the systematic study of 
many different CDW compounds^. We thus relate the 
C-u ~ contribution to strong pinning impurities22i, 

even though the microscopic nature of these impurities 
is not well understood experimentally (see Ref. |23 for a 
study of ESR spectroscopy in o-TaSa). 

The commensurate organic spin-Peierls com- 
pound (TMTTF)2PF6 showing slow relaxations^ 
contrasts with the inorganic spin-Peierls compound 
Cui_j;Zna;Ge03^2i2ii224^ showing antiferromagnetic 
ordering. We argue that the difference lies in the 
different nature of disorder. Substitutional disorder 
relevant to Cui-j^ZuajGeOa is qualitatively different 
from strong pinning impurities in CDWs and SDWs. 
The spin-Peierls compound Cui-xZnj^GeOa has a 
fast dynamics, with a "microscopic" time presumably 
comparable to the one of spin glasses (tq — 10~^^ sec) 
whereas in CDW and SDW compounds we have tq ~ f 
sec for the thermally activated process. This indicates 
that rather different mechanisms are at work, identified 
here are substitutional or strong pinning disorder. We 
generalize to the incommensurate case the model of 
substitutional disorder introduced in Ref. |3 In this 
model the solitons are due to domain walls between two 
degenerate ground states since the imp urity site can be 
removed from the chain (see Ref. |3J and section , 
therefore leaving randomly distributed domain walls 
in the chain from which the impurities sites have been 
removed. The specificity of this model (as opposed to the 
strong pinning model) is that there are no metastable 
states of bisolitons like in the strong pinning model. 
The substitutional impurities do not contribute to the 
slow dynamics of bisolitons in the strong pinning limit 
but are expected to contribute to collective pinning of 
the phase of the density wave. However we consider 
here temperatures much lower than the glass transition 



temperaturi 



.35.36 



so that interactions among solitons are 



the only remaining collective effects. For substitutional 
disorder we find interactions among solitons due to 
Friedel oscillations. We obtain similarly to Ref. a 
power-law specific heat C„(T) ^ T" and a susceptibility 
x(r) ^ T^-i+a agreement with existing experiments 
on the CDW o-TaSgSI. 

The article is organized as follows. In section [H] we 
investigate a classical model of collective effects in a dis- 



ordered CDW. Quantum effects are investigated in sec- 
tion HTH Final remarks are given in section HVI 



II. CLASSICAL LIMIT (STRONG PINNING 
IMPURITIES) 

A. Hamiltonian 

Let us start with a classical model of disordered 
CpVyi5,i6,i7,i8,i9,20,2i,22,23,24^ ^0 derive the ID projec- 
tion of the Hamiltonian of the phase of the CDW in the 
mean field approximation we follow the recent review by 
Brazovskii and Nattermann24 and consider a system of 
coupled chains with a phase ipniu) in chain n [y is the 
coordinate along the chain axis). 



An 



dy 



(2) 



dy [I - cos {(pn{y) - ifniiy))] (3) 



(«) 



1 — cos 



+ (4) 



where the sum in the last term runs over all impurities, 
vp is the Fermi velocity along the chain axis, Wm,n cor- 
responds to the commensurate energy or interchain cou- 
pling, V^"'^ is the pinning energy of the impurity number 
i in chain n, and Q = 2kF is the wave vector of the 
CDW. Assuming dilute impurities we suppose that the 
chain n = with (po{x) = ip{x) is coupled to neighboring 
chains with ipm{x) = 0. We arrive at the effective ID 
Hamiltonian 



n 



hvp 

An 



+ w / dy[l~ COS ipiy)] 



^V,[l-cos{Qy., + ^{y^))], 



(5) 



Metastable states due to the competition between the 
commensurate potential and the pinning energy were 
first discussed by Abe^^ in a different approach. Soli- 
tons and the transition to a 3D density wave glass were 
discussed by Fukuyama^^, and bisolitons were discussed 
by LarkinSS and Ovchiniko\i^. 



B. No impurity: 27r-solitons 

Without impurities there exist solutions minimizing 
the energy (jSJ in which the the phase winds by ±2n 
within a length 



(6) 
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where the soliton is centered at xq 
The width of the sohton is 



hvp 

2'KW 



zb^n (tan(V'/4)). 



(7) 



The energy Es of the sohton defined by Eq. 10 is equal 
to 4u;$i^: 



Es 



2hVFW 



(8) 



The energy Es can be viewed as the phase ordering 
temperature. Nevertheless thermodynamic equilibrium 
is not reached since the phase is frozen at a temperature 
larger than Es due to collective pinningSiiSi^. At the 
very low temperature considered here the relevant exci- 
tations are local deformations of the CDW in the form of 
bisolitons. 



C. One impurity: bisolitons 

Minimizing the energy with respect to ^{y) leads to 
the exact expression of the phase profile of the bisoliton 
associated to an impurity located at position yi with a 
pinning potential Vi. 



where ^ is such that 



tan I — ) exp 



\y-yi\ 



sm 



2hvp 



sin (ai + ip), 



(9) 



(10) 



with ai = Qyi- The bisoliton has a decay length given by 
Eq. 10. There can be several solutions to the matching 
equation (fTn|l . corresponding to metastable states sepa- 
rated by energy barriersSSiS^. 



D. A finite concentration of impurities 

1. Transition to a 3D density wave glass 

The problem of a finite concentration of impurities 
was treated in Ref. 22 within dynamical renormalization 
group (RG) describing a quench from high temperature. 
The first step was to consider two bisolitons at distance 
R. If R is much larger than ^ then the two bisolitons 
are independent from each other and have independent 
dynamics. If R is much smaller than ^ there is a single 
bisoliton pinned by two impurities, with a longer relax- 
ation time than two independent solitons. The case of 
intermediate values of R was treated approximatively by 
an exponential interpolation between the limiting cases 
i? <C ^ and R ^ In the following we consider a sim- 
plified model where the pinning energy is additive for 
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FIG. 1: Schematic representation of the spatial variation of 
the phase '^{y) along a chain, (a) A 2n soliton varying over a 
length ^. (b) A bisoliton— (superposition of a 27r-soliton and 
a 27r-anti soliton) varying over a length ^ generated by one 
strong pinning impurity, (c) Randomly distributed bisolitons 
generated by the clustering of strong pinning impurities for 
^ 1- (d) Density wave glass for a;^ > 1. The filled squares 
represent the impurities. 



all values of R smaller than ^, whereas the two bisolitons 
are independent from each other for all values of R larger 
than ^. Considering not only two impurities but a finite 
concentration of impurities distributed at random in ID, 
this defines clusters of impurities: two neighboring impu- 
rities at distance R belong to the same cluster if i? < ^, 
and belong to two different clusters if i? > f . 

The probability to find A^imp > 1 impurities in a given 
cluster follows the exponential distribution 



V{Ni^p) = exp {~xO [1 " exp {-x^] 



(11) 



where x is the impurity concentration. The average 
number of impurities within a given cluster is equal to 
exp (x^). For > 1 the number of impurities in a given 
cluster becomes very large so that the phase is effectively 
frozen because of large energy barriers. In this regime 
the system becomes a 3D density wave glass^^ and our 
treatment (valid if x^ < 1) based on the crude ID mean 
field model breaks down. 

We suppose in the following that the only 3D effect 
is this Larkin-Ovchiniko\*2SiS level splitting mechanism. 
We start from the limit of very dilute impurities and in- 
crease progressively the concentration. The crude mean 
field model is expected to describe well localized excita- 
tions in the dilute limit < 1 but at large scale and tem- 
peratures higher than the one considered here the system 
becomes a 3D elastic medium with different properties 
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that we do not discuss in the following. For instance the 
impurity perturbation decays as a power-law if I? > 1 
(see Ref.|23). 

2. Multi-impurity energy landscape 

We suppose that the phase is almost constant and 
equal to t/j in the middle of a bisoliton and note X = 
tan (■0/4) . The elastic energy corresponding to the sec- 
ond term of Eq. jSJ is independent on the number iVimp 
of impurities because it arises from the two sides of the 
bisoliton profile where dip{y)/dy is important. The pin- 
ning energy is additive because each impurity brings its 
own pinning energy. 

The energy landscape of a cluster of iVi,„p impurities 
generalizing Ref. is given by 



E(X) 







/ > 






/ ^ 


\e+av. 

\ / 


E 







E(X) 

A 



E -lA^ 





E 4AE 





(a) 



(b) 



FIG. 2: Schematic representation of the energy landscape in 
the commensurate case (a) and in the incommensurate case 
(b). The ground state is at energy Eq. The metastable state 
is at energy Eo + AE, where AE is the splitting. The unstable 
"bounce" state is at energy Eo + AV where AV is the energy 
barrier. 



E{X) = 16< 



1 + X^ 



Wimp |- ^ 



x^ 



x^ 



COS 



2X 



2 J 1 + X^ 



(12) 
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where the elastic and commensurate energies in the first 
term are much smaller in magnitude than the pinning 
energy, but play a role in lifting the degeneracy of the 
effective two level system. 

To illustrate Eq. I|12|) we show that minimizing with 
respect to X for impurities at a distance much smaller 
than ^ is equivalent to solving the sine-Gordon equation 
in the presence of the pinning term. Let us consider iVjmp 
impurities at positions j/i, yNi„-,p and denote by -01, 
ipNitn the value of the phase at the points j/i, yNi^^ ■ 
The solution with the appropriate boundary conditions 
is 



tan 



v(y) 



-01 \ ( y ~ yi 

tan ( — 1 exp I 



for y < j/i, and 

'f{y) 



tan 



tan [ — ) exp 



y - yNin 



(13) 



(14) 



for y > yNin-ip- The derivative of f{y) is discontinuous at 
the position of the impurities: 

+s dip _ 27rVfc . 

a^(%)-a^(y.) = 7— + (15) 

Assuming that tp is almost constant we obtain 

Wimp 



sm 



|) = ]S:E^^«-K + ^)' (16) 



hvp 



k=l 



generahzing Eq. l(TT)|l . Solving dE(X)/dX = with 
E{X) given by Eq. and X = tan (0^/4) leads directly 
to Eq. lfTC|) . 



E. Properties of the energy landscape 

The energy landscape E{X) given by ((T^ describes a 
ground state, separated from a metastable state by an 
energy barrier. There might be more than two energy 
minima for some realizations of disorder^^-^'^. In this case 
we restrict the energy landscape to the ground state and 
to the energy minimum separated from the ground state 
by the lowest energy barrier. We note /S.E the differ- 
ence between the energies of the metastable state and 
the ground state and /AV the energy barrier (see Fig.[2J). 



1. Commensurate case 

In the dimerized case {Q ~ a — Q) AE is equal 
to zero since the energy landscape is symmetric under 
a change of sign of the "coordinate" X. This can be 
seen from Eq. H12|l by noting that cos(ai/2) = ±1, 
sin(ai/2) = if Qxi — 27rn, with n an integer, and 
cos (ai/2) = 0, sin {ai/2) — ±1 if Qxi = {2n + l)7r, with 
n an integer. The degeneracy can also be seen directly 
from Eq. H16|) that, in the commensurate case Q = n, 
becomes 



,'0 

cosl- 



27r£ 



iVi, 



fe=l 



Vfc cos ak 



(17) 



If is a solution then — is also a solution, explaining 
the degeneracy. The degeneracy can be also understood 
by noting that the transformation ip{x) —ip{x) is a 
symmetry of the pure system, that is preserved by the 
pinning term for commensurate impurities but not for 
incommensurate impurities. 

In the commensurate case it is thus not possible to 
communicate energy over long time scales to the effective 
two-level system by increasing temperature. The classi- 
cal model predicts no slow relaxation at all in the com- 
mensurate case whereas in experiments there exists slow 
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Energy splitting distribution 
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FIG. 3: Distribution of level splitting /S.E for the model with 
random phases discussed in section Hi Dl The impurities are 
distributed at random on a chain of length 5. The impurity 
concentration is such that the average number of impurities 
is Ni-mp = x(^ = 1. The curves corresponding to 2;^ = 0.5 and 
a;^ = 5 are superposed on the curves corresponding to a;^ = 1. 
The pinning potentials Vi are uniformly distributed in the in- 
terval [0, Vo] with Vo indicated on the figure. The distribution 
of level splitting normalized to AEs is almost independent on 
and V(i, chosen such that AEs ^ Vq. The dashed line 
corresponds to ^(2) — l/[2-\/l — z\, with 2 — AE/AEs (see 

Eg. nil. 



relaxation, even though faster than in the incommensu- 
rate compounds^. Adding quantum tunneling between 
the two energy minima of the energy landscape can gen- 
erate two non degenerate energy levels corresponding to 
symmetric and antisymmetric wave functions, therefore 
restoring a finite heat response. 



2. Energy splitting distribution in the incommensurate case 

The distribution of energy splitting is shown on Fig. |31 
The distribution of splitting is close to 



Pi 



1 



AEs) 



(18) 



The most probable level spacing is close to AEs, inde- 
pendent on the value of the pinning potential and on 
the number of impurities involved in the bisoliton, show- 
ing that the shape of P{AE /AEg) is almost unchanged 
when the concentration of impurities increases. The en- 
ergy splitting distribution (|18|) can be understood in the 
case of a single soliton by noting that the pinning energy 
is much larger than the elastic energy so that the elastic 
term can be treated as a perturbation. There are four 
values of X minimizing the pinning energy: 



X, 



— sin {a/2) 
cos (a/2) 



(19) 



X 



e -I- cos (a/2) 
sin (a/2) ' 



(20) 



with e = ±1. The pinning energy of the solution H19|) 
is i5pin(^o'^'') = — 2y, and the pinning energy of the so- 



0. The elastic energy of the 



lution |2nj is Epin{x[''') 
solution X^'^ isE^iiX^'^) = 2Es{l-esiii{a/2)). The en- 
ergy sphtting i?ci(-Yo '') — Eci{Xq'^^) is then distributed 
according to Eq. (|18|l since a is uniformly distributed. 

The existence of the upper bound AEs in the energy 
splitting distribution is compatible with the experimen- 
tal observation of a high temperature tail ~ 1 /T^ of a 
Schottky anomaly in the equilibrium specific heat, with 
a well-defined level splitting. Experimentally the level 
splitting AEq of a two-level system is related to the tem- 
perature Tmax of the maximum of the Schottky anomaly 
by the relation AEq — 2.5/cBTinax- It was shown exper- 
imentally that Tinax < 30 mK^ so that A^o = iEs < 
100 mK. The existence of a well-defined level splitting in 
experiments is a universal property, valid for commensu- 
rate and incommensurate systems, and independent on 
the value of the CDW or SDW critical temperature that 
can vary by more than one order of magnitude from one 
compound to the other (Tpoieris = 218 K for o-TaSa, 
TsDW = 12 K for (TMTSF)2PF6, Tgp = 15 K for the 
spin-Peierls compound (TMTTF)2PF6, Tap = 13 K for 
the antiferromagnet (TMTTF)2Br). This is compatible 
with the fact that AEq is related only to the strength w 
of interchain interactions and the Fermi velocity vp: 



AEq = AEs = 



16 

s/2^ 



VpW, 



(21) 



where vp has the same magnitude for several com- 
pounds (for instance vp = 0.86 x 10'' cm. see" ^ in 
(TMTSF)2PFg^) and interchain couplings are also ex- 
pected to take similar orders of magnitude, from what 
we deduce that AEq is almost identical in all samples, 
independent on the charge or spin gap. The energy AEq 
is proportional to the elastic energy Eg of a 27r-soliton 
in the absence of impurities, not to be confused with the 
activation energy Ea of the order of 100 K probed in 
transport experiments^^. There is thus no inconsistency 
in the difference between the orders of magnitude of Ea 
and £^5. 



3. Distribution of energy barriers 

The distribution of energy barriers P{AV/Vo) is well 
fitted by an exponential function (see Fig. ^ . The aver- 
age relaxation time defined as 



To J d{AV)P{AV)T,exp{AV/T), 



(22) 



with T the temperature diverges at a finite temperature, 
like in the REM- like trap model^^. However this model, 
useful for discussing waiting time effects, over evaluates 
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FIG. 4: Distribution of the dimensionless energy barrier 
AV/Vo in the cluster model. We use — 0.5, 1, 1.5, 
kp = 7r/2 + 0.02. The pinning potentials Vi are uni- 
formly distributed in the interval [0, Vb] with Vo — 50. 
The solid lines correspond to the fits P{AV/Vo) = 1.45 x 
exp(-0.76 X AF/Fo) for x£_ = 0.5, P{AV/Vo) = 2.25 x 
exp(-1.04 X AF/Fo) for = 1, and P{AV/Vo) = 6 x 
exp(-1.61 X AV/Vo) for x^ = 1.5. 



FIG. 5: Variation of the out-of-equilibrium total specific 
heat as a function of log^Q(r/4i5s), for different values 
of tm/ro: logio(f„/ro) = -2 (□), logio(t„/ro) = -1.2 
(■), logio(U/ro) = -0.4 (o), logio(t™/ro) = 0.4 (.). 
logj^Q(t„/ro) = 1.2 (A). We used a:^ = 1. The pinning poten- 
tials Vi are uniformly distributed in the interval [0, Vb] with 
Vo/Es = 12.5. The specific heat for long t^/To follows a l/r^ 
behavior for T/AEs > 1, where Es is the energy of a soliton. 
The solid line represents the Cv ~ behavior. 



glassiness compared to the dynamical RG already dis- 
cussed in Ref. [ij (see section III Gp . released divided by the temperature variation: 

r<tot(, ^, U{tu,,t^,T + AT)~U{0,0,T) 
F. Out-of-equilibrium specific heat ^^^^^J /^j^ • ^'^^l 



To evaluate energy relaxation in the incommensurate 
case, we restrict the energy landscape to the ground state, 
the metastable state and the unstable "bounce" state at 
the top of the barrier. We are left with a two-state trap 
model with energies —Ei and —E2 (with Ei > i?2). Ei is 
equal to the barrier AV calculated in the preceding sec- 
tion and E2 is equal to the barrier AT^ minus the splitting 
AE. The unstable state is at zero energy. We note Pi{t) 
(P2(i)) the probability to be in state "1" ("2") at time 
t and note Pi,2(i) = Pi,2(0 exp (--Bi,2/2T). The evolu- 
tion of the probabilities is given by Glauber dynamical 



To 



dt 



Pi 
P2 



G 



Pi 
P2 



(23) 



with Gi.i = -expi-Ei/T), G2.2 - -exp(-£;2/T), 
Gi,2 = G2.1 = exp(-(i;i -f £:2)/2T). Experimentally 
the time Tq associated to bisoliton dynamics is of order 
of 1 sec for the thermally activated process. The time de- 
pendence of the occupation probabilities Pi {t) and P2 {t) 
is obtained by diagonalizing the 2x2 matrix G, from what 
we deduce the value of the energy U{tw,tyj + t,T) as 
a function of the waiting time t^, the time r elapsed 
since the waiting time and temperature T. The out-of- 
equilibrium total specific heat is defined as the total heat 



The dynamics equations can be solved exactly in the case 
of a single two-level system with energies —Ei and —E2 
and expanded in the high temperature regime T ^ \Ei — 
E2 1 , leading to the specific heat 



Gy 



{El - E2y 

1 — exp 



(25) 



To 



exp - — 



El 
T 



E2 

exp 1 - — 



proportional to l/T^, as expected for the high tempera- 
ture tail of a Schottky anomaly. Increasing the waiting 
time increases the energy transfered to the two-level sys- 
tem and therefore increases the amplitude of the 
term. 

The variations of the out-of-equilibrium total specific 
heat C^°* {tyj ,T) as a function of temperature T for dif- 
ferent values of the waiting time ty, are shown on Fig. [S] 
for a finite concentration of incommensurate impuri- 
ties. At temperatures larger than the maximal splitting 
AE ~ AEs the out-of-equilibrium specific heat follows 
a behavior. The out-of-equilibrium specific heat is 
strongly reduced as the waiting time decreases but still 
follows a behavior, in a qualitative agreement with 
experiments (see Fig. 3 in Ref. 0) ■ 
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FIG. 6: Spectrum of relaxation times deduced from energy re- 
laxation (see Eq.lEll). We use Vo/AEs = 12.5, T/AEs = 400„ 
and a;^ = 1. The pinning potentials Vi are uniformly dis- 
tributed in the interval [0, Vb]. The waiting times correspond 
to logio(t»/ro) = 2.4 (□), 7.2 (■), 12 (.), 16.8 (o), 21.6 (A). 
The solid line is a fit to logj^Q(P(logjg r)) = a + & logj^g '''i with 
h = -0.083. 



G. Spectrum of relaxation times 

Let U{tw,tw + T, T) be the energy of the two-level sys- 
tem at time t — t^+r, with a heat pulse applied between 
t = and t — tw The spectrum of relaxation times is 
deduced from U {t^ , t^j -f r, T) by assuming that 

U{t^,t^ + T,T)= J Pt^{\nT')exp{~T/T')d\nT'. (26) 

Replacing the exponential by a step-function leads to the 
relaxation time spectrum^ 

(lnr)oc ^^^^-^+"'^^ (27) 
omr 

We have shown on Fig. |^ the spectra of relaxation times 
Pt^(lnr) for — 1. Similar results are obtained for 
= 0.5 and — 1.5, with however a different value 
of the exponent of the power-law. The long time tail of 
the spectrum is a power-law, in agreement with experi- 
ments (see Ref . Il2|) . As the waiting time increases the 
power-law regime extends to ever longer times without 
limits, compatible with the existence of a genuine glass 
transition that is an artifact of this clustering model. 

For comparison we have also calculated the relaxation 
time spectra obtained from dynamical RG (see Fig. [7|). 
The calculation is identical to Ref. we suppose a 
quench from high temperature at time t — 0. The 
small time degrees of freedom are progressively elimi- 
nated, making bigger object that relax more slowly. In 
this approach the energy landscape of two impurities at 
distance R is interpolated between the two limiting cases 



10" 

10 15 20 25 
logio(x/xo) 

FIG. 7: Spectrum of relaxation times deduced from dy- 
namical renormalization group. We use Vo/AEs = 12.5, 
T/AEs = 400,, and = 1.5 (■), a;^ = 1 (•) and = 0.5 
(a). The pinning potentials Vi are uniformly distributed in 
the interval [0, Vb]. The solid line is a fit to logjQ(P(logjQ r)) = 
a+felog^o r, with b = -0.23 for = 1.5, b = -0.21 for x^ ^ 1 
and b = —0.17 for x^ — 0.5. We carried out simulations with 
larger system sizes and found practically no finite size effects 
on the data on this figure. 

i? <C ^ and i? ^ ^ so that the phase fluctuations are more 
important compared to the previous clustering model. 
Impurities can thus be depinned in sequence, not nec- 
essarily all together. We obtain also a power-law relax- 
ation but (i) the exponent of the power-law spectrum 
of relation times is different from the clustering model 
(smaller relaxation times are favored), and (ii) ageing 
is "interrupted" for the clustering model (there exists a 
maximal relaxation time). As expected the clustering 
model discussed previously over evaluates glassiness but 
is nevertheless useful for addressing qualitatively waiting 
time effects that could not be discussed within dynamical 
RG. Experiments favor interrupted ageing with respect 
to a genuine dynamical glass transition (see Ref. [13 for 
an analysis of experiments based on a REM-like model 
that is justified microscopically by the dynamical RG ap- 
proach) . 



III. QUANTUM LIMIT (SUBSTITUTIONAL 
AND STRONG PINNING IMPURITIES) 

A. Preliminaries 

We base the discussion in this section on the electronic 
part of the Peierls Hamiltonian. We discuss a model of 
disordered CDW there is only a charge sector, or a model 
of disordered SDW where there is only a spin sector. We 
expect this simplified models to capture part of the prop- 
erties of real compound in which there is a gap in both 
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a' 



(d) 



(e) 



p b P' 



p P' 







FIG. 8: Schematic representation of a dimerized CDW chains 
with open boundary conditions with one impurity (a) or two 
impurities (b). The double lines represent the strong bonds t+ 
e. (b) corresponds to R = yp —ya odd. We represented on (c) 
two substitutional impurities (black square) leaving unpaired 
fermions at sites a and /3. (d) corresponds to the chain where 
the impurity sites have been removed—, (e) shows the spatial 
variation of the phase corresponding to the chain (d). 



the spin and charge channels. The analogy between the 
CDW case with only a charge channel and the SDW case 
with only a spin channel is based on the fact that one 
dimensional tight binding models can be mapped onto 
XX magnets in one dimension through a Jordan- Wigner 
transformation. It was shown explicitely in the dimer- 
ized case that the properties of substitutional disorder 
are almost identical for the XX and Heisenberg models^ 
in spite of interactions between Jordan- Wigner fermions 
in the Heisenberg case. We expect this similarity to be 
valid also in the incommensurate case. 

Quantum effect were already discussed in Ref. |^ but 
here we work within a different microscopic model that 
was proved to be useful in the case of doped spin-Peierls 
system^:, and generalize it to the case of incommen- 
surate CDWs. Substitutional disorder in a spin-Peierls 
system (relevant to Cui-x'Z'^xGeO^Si^i'^^^) is qualita- 
tively equivalent to impurities that break the chains into 
finite length segments (see Ref. Is^ and section lTlI Cll . We 
have shown on Fig. |S1 a schematic representation of sub- 
stitutional impurities in dimerized chains. The phase of 
the dimerization changes by tt in the chain where the 
site of the impurity has been removed. For a CDW or 
SDW with a wave- vector Q = 2kp the phase changes by 
QflQ once the impurity sites have been removed, where 
ao is the lattice parameter. The absolute value of the 
amplitude of the CDW is supposed to remain constant. 
The solitons generated by the substitutional disorder that 
we consider here are thus phase solitons, connecting two 
degenerate ground states with different phases. One 
can cross-over continously from edge states in the lim- 
its ta^a' = tf3^i3' = to solitons Connecting two ground 



states as the parameters ta^a' and t/^^p' are increased. 

There is a flat classical energy landscape in the case of 
open chains for the model defined by Eq.jS)) in the limit 
w — 0. The classical Hamiltonian is just |(SJ| on an open 
chain, without the pinning term: 



n 



An 



dy 



dy J 



(28) 



w I dy[l- cos(^(?/)] 

'0 



that has no bisoliton metastable states like in the preced- 
ing section (the ground state is (p{y) = 0). The bound 
state associated to a single substitutional impurity in the 
quantum model defined by the Peierls Hamiltonian 



H = ^ [< -I- e cos {2kFyi)] [ct+ic, + c+q+i] 



(29) 



is thus at the lowest possible energy (exactly in the mid- 
dle of the gap). The variable t in Eq. H29() is equal to 
the average hopping amplitude and e is the amplitude of 
the modulation. The variable yi in Ea. H29|l is the coor- 
dinate of the site number i: yi = iaQ, with ao the lattice 
parameter. The bound state due to substitutional disor- 
der can be occupied by an electron and a hole with an 
equal probability, and there is a fast dynamics of the oc- 
cupation probability since there are no energy barriers in 
the classical limit. Like in spin glasses mostly governed 
by RKKY interactions it is expected that the "micro- 
scopic time" tq is extremely small (tq ~ 10~^^ sec in 
spin glasses). By contrast there is slow dynamics with 
To ~ 1 sec in the case of the organic spin-Peierls com- 
pound (TMTTF)2PF6, well described by the classical 
model discussed previously. 

The width of a soliton in the quantum limit is equal to 



2hvp 



x/A2 - ' 



(30) 



where E is the energy, A is equal to the Peierls gap 
and vp is the Fermi velocity in the absence of lattice 
distortions, ^o is different from the soliton or bisoliton 
width in the classical limit ^ ~ ^fivp jlnw discussed 
in the preceding section. In the classical model one 
should incorporate interchain interactions for the soli- 
ton or bisoliton to have a finite width whereas in the 
quantum model the zero energy soliton width is finite 
in the absence of interchain interactions. Experimen- 
tally the slow relaxation properties associated to the 
l/T^ tail of the specific heat are independent on the 
value of the gap A and of the transition temperature 
that can vary by more than one order of magnitude. In 
the case of the SDW (TMTSF)2PF6 the BCS-like re- 
lation A ~ l.TTc is well verified (T^ = 11.5-^12 K 
and A = 20 K from NMR,^). The Fermi velocity^^ is 
vp ~ 0.86 X lO'^cm.sec"-^ so that ^o = 2hvp/l\ ~ 90ao, 
with the lattice parameter ao = 7.3A. In the case of 
o-TaSa we do not have a precise data for vp and we 
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take vf = lO^cm.sec "'^ as a typical value. We have 
j.^(CDW) ^ 215K and Acdw = 3.73Tc ~ 780K so that 
^(CDW) _ -virith qq = 3.34A. However it is likely that 
there is also a smaller spin gap since the spin suscepti- 
bility follows a power-law. An energy scale E = 30 K 
can be estimated from the deviations of the power-law. 
Assuming that the power-law susceptibility is due to ran- 
domly distributed magnetic moments we identify the en- 
ergy scale E to the average exchange energy J^v — Aa:^ 
with = 0.3 estimated from the power-law suscepti- 
bility (see section lllTTT2|) so that ^^'^'"^ = 150-=- 450ao. 
^^"P'"^ in the SDW compound (TMTSF)2PF6 and o-TaSg 
are thus one order of magnitude larger than ^0 in the 
spin-Peierls compound CuGeOs. 



kp = 71/2 kp = jt/2 + 0.02 




01234 01234 

FIG. 9: Evolution of the bound state energy levels as a func- 
tion of the distance between the two pinning centers, in the 
commensurate case kp = 7r/2 (a) and in the incommensurate 
case kp = 7r/2 -|- Skp, with Skp = 0.02 (b). The energy is in 
units of A = \/2e and the distance is in units of ^0 ~ 2vf/ A. 
We used e/t = 10"^ and V/t = 10. 



B. Strong pinning impurities 



C. Substitutional disorder 



1. Dyson matrix 

The quantum model given by Eq. (|29|l can also be used 
to treat strong pinning impurities, not only substitutional 
disorder. The Hamiltonian is equal to H29I) . plus a term 
describing the impurity potential: 



(31) 



where the strong pinning impurities are located at ran- 
dom positions {j/i, ...,2/Ni„p}- We suppose that Vy^ = V 
is the same for all impurities and that V > Q. We note 
G{E) the Green's function at energy E in the presence 
of the pinning potential and g{E) the Green's function in 
the absence of pinning potential. The Green's function 
G{E) is obtained from inverting the Dyson matrix: 

E [-^-^^ + 9y..y. {E)Vy.] Gy,,y, (E) = 5,,,,^. (E), (32) 
fc=i 

where the Green's function gy^^y^iE) is given in Ap- 
pendix The bound state energies correspond to the 
poles of G{E). 



2. Two pinning centers 

In the commensurate case kpcio — n/2 the energy lev- 
els evolve smoothly as a function of the distance between 
the two impurities (see Fig. EKa))- In the incommen- 
surate case kptto = n/2 + Skp (see Fig. EKb)) the en- 
ergy levels fluctuate strongly as the distance between the 
two impurities is reduced. This shows that quantum me- 
chanical interactions among impurities enhance disorder 
effects in the incommensurate case but not in the com- 
mensurate case. 



1. Open chains in the dimerized limit 

Let us consider two semi-infinite dimerized chains end- 
ing at sites "a" and "a" (see Fig. |HKa)). We note 
to — ta.a the value of the hopping between sites "a" and 
"a" (equal to i -I- e or t — e) in the infinite chain. The 
Dyson equation relates the Green's functions gij of the 
infinite chain to the Green's functions Gij of the semi 
infinite chain: 



gaAE) = Ga.aiE)+tlGa,a{E)gaAE) (33) 
(E) = Ga.AE)+tlGa,a{E)gaAE), (34) 



with ga.a{E) = ga.a{E) 

tion of 1(221) and is 



S/(2tV2e2 - E^). The solu- 



Ga.a[E) 



t V2e2 - E^ 



/2 



E 



-l + r]a\ 1 + 



t^ 

'•0 



E^ 



t2 2e2 - i;2 



(35) 

with r^a = ±1. In the case of two weak bonds in the chain 
(see Fig.lHKb)) we find Ga,i3 = ga.p/V, with 

V = [i + tl ,^Ga,agc.,c.] [1 + tlpGh^g^A (36) 

~ ta,atb,l3Ga.,aGb,bga,pgp,a- 

If the separation between the boundaries is large enough 
the energy of the bound states can be expanded in the 
small parameter exp(— i?/^o) where ^o/o-o = 2^/2t/e is 
the zero energy coherence length. Since the bound states 
are close to the Fermi energy we also expand Ga.p in 
powers of If R — yf} — ija is even we find that one 
bound state is generated at one end of the chain. In the 
case R odd on Fig. |Hl-(b) we find two bound states at 
E'-'^^ — ±(e/3) exp (— i?/^o) corresponding to rja = rjb ^ 
— 1. In the ground state E^~^ is occupied with an electron 
and is occupied with a hole. The lowest neutral 

excitation corresponds to an electron in the level E^~^^ 
and a hole in the level E^^\ This is in agreement with a 
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R/i;o R/5o R";o 



FIG. 10: Evolution of the bound state energies as a function 
of the distance R = \ya — yb\ between the two impurities. 
We use kp = ir/2 + Skp, with Skp = (a), Skp = 0.02 (b) 
and Skp = 0.04 (c). We use the parameters e/t = 10~^. 
ia.ci'/^ = ^li.P'/'t = 1/2. The dashed hue is a fit to E/A = 
±6.26 cos {{5kF)R) exp {-R/^o)- 

previous model developed for the spin-Peierls compound 
Cui-jjZn^Oj"^^ as well as in a qualitative agreement with 
numerical simulations^ii^. 



2. Substitutional disorder in CDWs and SDWs 

We consider now a more realistic model of substitu- 
tional disorder, both for commensurate and incommen- 
surate CDWs. We suppose that two substitutional impu- 
rities "remove" the sites "a" and "b" from the chain and 
that second neighbor interactions ta.a' and t/3,/3' couple 
the right and left neighbors (see Fig.|SKc)). The model is 
similar to Ref . [13 proposed for Zn impurities in CuGeOa 
except that we consider here the incommensurate case. 
We calculate the poles and residues as a function of the 
distance R — \ya — yb\ between the two impurities (see 

Fig.Hni). We use kp = k^p^ + Skp, with = 7r/2 and 
show on Fig.^lthe Friedel oscillations of the bound state 
levels for odd values of R. The phase of the Friedel oscil- 
lations is shifted by tt/2 for even values of R. We deduce 
from the level crossings that there exits a change of sign 
in the hopping amplitude between the two solitons as R 
is increased. 

Similarly to a model of doped spin-Peierls system.'^^ 
substitutional disorder in CDWs and SDWs provides an 
explanation to the power-law specific heat ~ T" ob- 
served in disordered CDWs and SDWs, and to the sus- 
ceptibility x{T) ~ r"i+" observed in the CDW o-TaSj^^ 
(with the same value of a in the specific heat and sus- 
ceptibility experiments). In experiments on the SDW 
(TMTSF)2PFgiS, on changing the time constant in a 
heat pulse experiment by a factor of 100, there is an 
increase by a factor of 7 of the amplitude of the 
contribution whereas the contribution changes only 
by 20%, and not in a systematic manner. In the com- 
mensurate (TMTTF)2Br compound and also in a heat 
pulse experiment, an increase of the time constant by 
a factor of 3 results in an increase of the term by 
a factor of 3 and only an increase of the T" term by 
20%. This shows that the T" contribution to the spe- 



cific heat can be interpreted as an equilibrium property, 
which is the case if frustration due to Friedel oscillations 
can be neglected as in a dilute one dimensional model. 
Focusing on the SDW case we consider only exchanges 
J{R) between nearest neighbor solitons and disregard the 
changes of sign in J{R). The exchange distribution scales 
like P(|J|) ~ \J\~'^^^^°, where x is the concentration of 
substitutional impurities^. The energy of "active" pairs 
of spins is approximately given by f7 = /g JP{J)dJ, 
where T is the temperature so that the temperature de- 
pendence of the specific heat is given by ~ T^io ^ xhe 
susceptibility is approximately given by a Curie contri- 
bution for the fraction of "active" spin having an en- 
ergy smaller than T so that the susceptibility behaves 
like x{T) ^ T~^^'^^°, leading to a = x£,q. Experimen- 
tally a = 0.3 ±1.2 as mentioned in the Introduction. 
In the case of o-TaSs we estimate — 150 ± 450ao in 
the spin sector as mentioned previously and a — 0.3 ob- 
tained from specific heat and susceptibility. We deduce 
the concentration of intrinsic substitutional impurities 
Xint — 0.07 ± 0.2%. The upper bound is comparable 
to the nominal concentration of extrinsic Nb impurities 
Xcxt — 0.5% which might explain why the exponent a in 
the experiment is the same in the presence or absence of 
extrinsic impurities. 



IV. CONCLUSIONS 

To conclude we have investigated several factors in- 
volved in energy relaxation in disordered CDWs and 
SDWs. A first factor is the role of commensuration. In 
the commensurate case we find that the energy landscape 
is symmetric since the symmetry ip{x) —ip{x) of the 
CDW Hamiltonian is preserved by the impurity poten- 
tial. As a consequence the two energy minima are degen- 
erate, a property that can be obtained by a direct solu- 
tion of the sine-Gordon equation. We have shown that 
the degeneracy exists for an arbitrary number of impu- 
rities contributing to pinning the bisoliton. Experimen- 
tally there exists also slow relaxation in commensurate 
systems, even though faster than in the incommensurate 
case. The fully classical model might be too schematic 
since the exact degeneracy of the classical model is lifted 
by quantum tunneling. Nevertheless it succeeds to ex- 
plain the Cy ~ T^^ specific heat and the waiting time 
dependence of the prefactor in the incommensurate case. 
We have proposed a model of clustering in which the 
pinning energy is additive if two impurities are at a dis- 
tance smaller than the width ^ of the soliton. The bar- 
rier distribution is exponential, with therefore a power- 
law relaxation compatible with experiments (see Ref. il 21') . 
However bisolitons can be depinned in sequence, not si- 
multaneously. This is described qualitatively by the dy- 
namical renormalization groupie. Within this treatment 
we obtain also a power-law relaxation but with an upper 
cut-off in the spectrum of relaxation times, in agreement 
with experiments. The clustering model is nevertheless 
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useful for addressing qualitatively the waiting time effect 
in the specific ageing protocol used in experiments. 



APPENDIX A: GREEN'S FUNCTIONS OF A 
CHARGE DENSITY WAVE 



Finally we discussed another important factor: the na- 
ture of disorder at the level of a single soliton and the 
possibility of quantum interactions among solitons that 
might explain the differences between the spin-Peierls 
compound Cui-^jZuj-GeOa that does not show slow re- 
laxation and orders antiferromagnetically, and the spin- 
Peierls compound (TMTTF)2PF6 that shows slow relax- 
ation without antiferromagnetic order. For this purpose 
we have generalized to incommensurate systems a model 
of substitutional disorder originally introduced for spin- 
Peierls systems^. We first treated the case of a strong 
pinning impurity potential and found that the energy 
spectrum of a single quantum mechanical soliton is con- 
stant as the position of the soliton varies along the chain 
in the commensurate case, whereas it varies in the incom- 
mensurate case. Interactions among solitons enhance dis- 
order effects in the incommensurate case. We thus find 
a qualitative difference between the commensurate and 
incommensurate cases in the quantum limit. For substi- 
tutional disorder in the quantum limit we find interac- 
tions among solitons due to Friedel oscillations. This can 
explain the experimentally observed power-law contribu- 
tion to the specific heat in inorganic CDWs and SDWs 
in the limit of weak interchain couplings. 

The final picture for organic SDWs and CDWs is a co- 
existence between strong pinning and substitutional dis- 
order as well as a a coexistence between classical and 
quantum effects. Further investigations would require a 
model of collective effects interpolating between the clas- 
sical and quantum limits. Another ingredient that we 
did not discuss in detail is metallic island formed around 
impurities^ that would lead to a phenomenology close 
to that of substitutional disorder because both generate 
states at the Fermi level for an isolated impurity, that 
can interact through Friedel oscillations. 

Finally, we have developed here the point of view of 
including collective effects from the strong pinning limit. 
The specific heat including quantum effects in the weak 
pinning regime was discussed recentlji^. 



1. Canonical transformation 

We use the Peierls Hamiltonian for charge degrees of 
freedom, the electronic part of which is given by Eq. H29|l , 
where spinless fermions jump between neighboring sites 
on a ID chain. We suppose that the lattice is frozen and 
therefore we do not include in the Hamiltonian the term 
due to the lattice deformations. The Hamiltonian (|29|l is 
diagonal in terms of the operators ^k,R and ^k,L, where 
R and L label right and left-moving fermions: 



n 



E 

k 



[Ek., 



RlkRlk,R + EkXlkLlk,L 



with 



-k.R 



^k,L 



'^k ^k+2kF,R 



and 



EkM = 2tcos(fcao) 



1 



(Al) 

(A2) 
(A3) 

.(A4) 



At sin (kpao) {k — kF)ao 

^ 1 

Ek,L = 2tcos(fcao)+ . \ 11.^1. A ^A5) 

At sm (Ki?ao) (fc + kFjao 

The coefficients i3^^^ and b'^'' are given by 



£ exp(zfcfao) 

At {k- kF)ao 



B 



(L) _ e exp {-ikpao) 



At [k + fcF)ao 



(A6) 
(A7) 



The normalization coefficients are given by J^j/^' 



l/y 1 + jiB^''^''^'' p. The Green's functions deduced from 
the spectral representations are given in Appendix El 

2. Green's functions 
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The advanced Green's function defined as 

gyuvA^iM) = -iO{ti - t2){{cy^{ti),Cy^{t2)}) de- 
composes into the sum of the four right (R) and left (L) 
combinations: 

9y^,yAtlM) - 9v;''y,{tiM) + 9^;^y,{h,t2) (A8) 
where the "RR" Green's function is defined by 

<l(ii'^2)= E e''''y^e-^^^y^{{ckMti)^ck,Ah)}). 

ki,k2 

(A9) 
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and similar expressions are obtain for the "RL", "LR" 
and "LL" Green's functions. The spectral representation 
of is given by 



1 



\B 



(fl)|2 



E — Ekji E — Ek-2kFX 



(AlO) 
,(A11) 



and similar expressions are obtained for the three other 
Green's functions. After performing the integral over 
wave vector in the spectral representations we obtain 



1 



V2€ 



sm < Lp 



E 



1 



1 



At sin k p 
\/2e 



3 2tsinkF V2e2 - 



R\e^v{-R/S,{E)) 



(A12) 



3 2tsinfcF V2e2 - E^ 



sm < tp 



E 



At sin (fci?ao) 



R\exp{-R/^{E)), 



where R = 2/i ~ 2/2 is positive and ^{E) ~ 
At sin {kpao) /V'^e'^ — E"^ is the coherence length at a fi- 
nite frequency. The sum of the "RL" and "LR" Green's 
functions is given by 



V2 1 V2e 

~ 2t sin {kpaa) ^/W^ 



(A13) 



exp(-i?/e(i?)) 



fcF + 2^ - 



kp 



E 



At sin (fcFOo) 



R - ikpyi 



cos 



kp 



E 



At sin (kpao) 



R-2kpy 



We obtain similar expressions for R < 0. 
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